Wide-band detection of the third moment of shot noise by a hysteretic Josephson 

junction 
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We use a hysteretic Josephson junction as an on-chip detector of the third moment of shot noise 
of a tunnel junction. The detectable bandwidth is determined by the plasma frequency of the 
detector, which is about 50 GHz in the present experiment. The third moment of shot noise results 
in a measurable change of the switching rate when reversing polarity of the current through the 
noise source. We analyze the observed asymmetry assuming adiabatic response of the detector. 



Studies of shot noise in mesoscopic conductors are 
presently of great interest, revealing detailed informa- 
tion on microscopic mechanisms of electronic transport 
0, S Q- The theory of full counting statistics (FCS) 
of electrons 0] determines the probability distribution 
of current fluctuations and its nth order moments. Ex- 
perimentally, however, measurements of higher moments 
beyond the variance remain difficult especially in the in- 
teresting high frequency regime, mainly because of weak 
signals and demanding filtering requirements. The first 
experiments on the third moment of shot noise of a tunnel 
junction were performed on samples connected remotely 
to the detector [B|, @] . Yet the most natural way to in- 
vestigate noise in nano-structures would seem to be via 
an on-chip detector. In recent measurements of Gustavs- 
son et al. [7( and of Fujisawa et al. [8|], single electrons 
tunneling through quantum dots could be counted di- 
rectly. The noise correlations were observed up to the 
third order over a bandwidth of a few tens of kHz. The 
latest remarkable results on observation of the fourth and 
the fifth moments of current through a quantum dot [9( 
demonstrated on-chip detection of FCS in the low fre- 
quency regime. As regards to high-frequency on-chip de- 
tectors of non-Gaussian noise, a Josephson junction (JJ) 
appears to be a very attractive one according to various 
theoretical predictions [ill . Hil Il3j |. On the experi- 
mental side, the influence of non-Gaussian noise on the 
conductance of a Coulomb blockaded JJ has been mea- 
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FIG. 1: Measurement scheme for on-chip detection of noise. 
A Josephson junction detects the current fluctuations aris- 
ing from the noise source. Only heavily filtered DC lines 
(through cold resistors, line inductances and thermocoax ca- 
bles) connect the setup with the room temperature amplifiers 
and current sources. On top right we show IV characteristics 
of the detector junction taken with different values of average 
current In through the NIS noise source. The IV of this NIS 
source is depicted in the lower right graph at small currents. 
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sured in [14] . In this Letter we present an experimental 
observation and a theoretical interpretation of the influ- 
ence of wide-band third order fluctuations on escape of 
a hysteretic JJ threshold detector from its zero voltage 
state. 

Our analysis of the experiment is based on the follow- 
ing model. Let us assume for a moment that the Joseph- 
son potential of the detector is strictly harmonic with os- 
cillation (plasma) frequency uj p — y^^f [1 - {I / Ic) 2 } 1 ^ 4 , 
where Iq and C are the critical current and capacitance 
of the junction at bias current /, respectively. Then, 
on the level of the second moment, noise affects the es- 
cape of a JJ predominantly by promoting resonant ac- 
tivation at oj p such that shot noise induced transitions 
between adjacent energy levels j and j — 1 in the po- 



tential are given by = 2h ^ c Si(—ll> p j 

and — 2 hL c '-M+^p) (relaxation), where S/(u>) 

is the spectral density of noise at (angular) frequency ui. 
These rates lead to dynamics in the potential with effec- 
tive temperature [l5| 



T* ~ /^p/[2/c B arcoth(l 



QF\h 
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Here we have assumed that the bath temperature is 
T <C hujp/k-Q, i.e., the detector is in macroscopic quan- 
tum tunneling (MQT) regime in the absence of shot 
noise, and that the voltage V across the scatterer is 
large enough such that we are in pure shot noise regime: 
|eV| > k B T,hu! p . In Eq. (Q} Q is the quality factor 
of the junction at plasma frequency, and /n is the aver- 
age value of the noise current In through the scatterer 
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junction with Fano factor F 1 1 61) . Unlike in the earlier 
experiments by some of us [1 51 ] . in the present measure- 
ments we can determine Q independently by detecting 
the crossover from escape dynamics to under-damped 
phase diffusion (UPD) [l7| as described in the experi- 
mental part below. This leaves no free parameters to 
determine the consistency between the measured escape 
rate 



TABLE I: Parameters of the samples. 



r 



-^exp(-A£//fcBT*) 

Z7T 



(2) 



and T* given by independent measurements of the pa- 
rameters in Eq. Q. In © AU ~ *&Ej(l - I/I c f /2 
is the barrier height of the potential well with Josephson 
energy E 3 = M c /2e. 

The linear coupling to current fluctuations SI in the 
Josephson potential is of the form —^-8I(p, where <p 
is the phase operator. It is proportional to the sum 
of the creation and annihilation operators of the har- 
monic oscillator, and therefore, in contrast to the influ- 
ence of the second moment, the third order fluctuations 
do not induce transitions between adjacent levels in a 
harmonic well 
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1181 ] . Therefore, up to the third or- 
der, the influence of resonant activation in a harmonic 
potential is exclusively that due to the second order 
fluctuations. There are weak corrections to the results 
above, if we allow the potential to be anharmonic. Yet 
in the case of the third order fluctuations such correc- 
tions vanish unless the third order spectral densities are 

" 3. There- 
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assumed to have frequency dependence 
fore we assume in what follows that the third order ef- 
fects come from non-resonant fluctuations only, at sub- 
plasma frequencies. Throughout we ignore the effects 
related to higher than third order fluctuations. Under 
these conditions, in analogy to adiabatic Gaussian noise 
[li|, the contribution of the third moment to the escape 
rate from the well is determined by the instantaneous 
value of the fluctuating bias current. The asymmetry 
Ar/r avc = «T+) - (T-))/[«r+) + (T-))/2] between es- 
cape rates (r 1 * 1 ), averaged over the adiabatic fluctuations, 
at different polarities of either pulse or noise currents can 
then be written as Ar/r avo ~ -|(M) 3 (<Jj3). Here 
B is the exponent in the expression of the tunneling rate 
and (61^) is the third moment of current fluctuations at 
the detector. For thermal activation of Eq. |J2J), with 
B = AU/kBT*, one then obtains 



Ar/r a 



W2 h 3 
3 K 2e> y 



I/Ic) 



3/2 



(k B T*) 3 



(3) 



u! p is the bandwidth 



Here, (61*) = (^) 2 S 3 , where Au 
of adiabatic fluctuations, and S3 = F 3 e 2 I^s is the low 
frequency limit of the third order spectral density with 
Fano factor F3. 

The measurement scheme is shown in Fig. [T] The de- 
tector is an Al/AlOx/Al J J with an area of ~ 2 /mi 2 



Sample 


detector 
Ic (M) G (fF) w p /27r (GHz) 


noise junction 
i? N (kfi) C N (fF) 


NIS 
SIS 
REF 


2.88 80 53 
2.51 80 49 
2.36 75 49 


15.4 8 
1.9 6 




and it stays initially in a superconducting state. Another 
tunnel junction, the noise source, is biased with a current 
±Jnj such that it is driven far from equilibrium into the 
shot noise regime. All samples were fabricated by elec- 
tron beam lithography and shadow evaporation and they 
were measured via filtered signal lines at a bath tempera- 
ture of ~ 30 mK. In order to detect only the fluctuations 
of In, the balancing current =p/fj of opposite polarity is 
applied such that no DC-component due to In passes 
through the detector. Owing to the small magnitude of 
the third order fluctuations, careful DC-balancing is ex- 
tremely important. The residual non-balanced current 
across the detector is measured repeatedly by a low in- 
put impedance current amplifier. However, the resulting 
correction never exceeded 1 nA during the measurements. 
The on-chip inductors in Fig. [T] are long (few mm) and 
narrow (2 /im) superconducting lines. In all the sam- 
ples, the detector was connected to a large contact pad 
on the chip (few mm 2 in area), which served as a capac- 
itive short to ground at high frequencies. To probe the 
fluctuations of 7pj, we applied trapezoidal current pulses 
of height I through the detector with At = 100 fj,s - 1 
ms duration, and with a 1 ms delay between two pulses. 
Typically 10 3 — 10 4 pulses at each value of I were re- 
peated, and the escape probability P was obtained as the 
fraction of pulses leading to an escape from the supercur- 
rent state. When increasing the average value of In , shot 
noise enhances the escape rate, leading to suppression of 
the escape threshold current [15] . To extract information 
about the third moment, we applied four different com- 
binations of current polarities, which we call / ++ , / , 
1^ , and / *~ , where the superscripts refer to the signs of 
I and In, respectively. The first two combinations should 
lead to identical escape characteristics, which should dif- 
fer, in the presence of odd moments, from the last two 
identical combinations. We normally measure the sepa- 
ration of the histograms, I + — I~ , corresponding to equal 
escape probability P in the two pulse/noise current direc- 
tions. By a simple geometric argument, assuming linear 
dependence between escape probability and current over 
a short interval (weak third order effects), we obtain 

31 

i + - 1-^(1- P) ln(l - p)(^)Ar/r avc . (4) 

We present data of three samples, see Table I. Sam- 
ple NIS had a normal mctal-insulator-superconductor 
Cu/AlOx/Al noise source. In Sample SIS the noise 
source was another Al/AlOx/Al J J, but with a much 



3 




1 ' I 

:(C) 









0.51 



0.54 



0.57 



0.87 0.88 0.89 
I/L 



FIG. 2: Escape characteristics of Sample NIS. (a) Histograms 
taken with the four different polarity configurations of noise 
and pulse currents, and at several values of noise current: 
|/n| = 0...3.8 fiA from right to left, (b) Escape rates cor- 
responding to |/n| = 2.2,2.4, and 2.6 fiK from right to left 
displaying the splitting of data with I ++ and I~~ currents 
from those with J" 1 and I~ + currents, respectively, (c) Es- 
cape rates with In = in MQT regime. 



smaller critical current than in the detector. In the refer- 
ence sample, Sample REF, the noise source was replaced 
by a superconducting line. In Table I denotes the 
normal state resistance of the noise junction, and w p is 
taken at I = 0. 

Here we describe data of Sample NIS unless otherwise 
specified. Figure [T] shows current- voltage characteristics 
of the detector and the noise source junction. The detec- 
tor IVs have been measured at different values of noise 
current. Besides the standard hysteretic character and 
suppression of switching threshold upon increasing 7n, 
there are two important features to observe here, (i) 
The IVs are vertical in the supercurrent branch, such 
that V < 0.2 /iV up to the switching current. This 
means that the rate of phase diffusion events from one 
well to another is / 2Tr = 2eV/h < 100 MHz. (ii) The 
gap voltage is the same at all values of noise current up 
to 5 /iA; this is higher than all those currents used for 
further analysis of data. Thus there is no excess over- 
heating of the bath and the only relevant temperature 
is then T* . Figure [2] shows escape characteristics with 
rates (T^) that are related to the measured escape prob- 
abilities P ± via P ± = 1 — e~^ r ^ At . Here the current 
pulses were At = 480 /is long. The measurements were 
performed at various values of In using the four polarity 
configurations. Suppression of the mean bias current and 
change in the width with increasing In is evident from 
the histograms in Fig[2](a). Moreover, as to the different 
polarity configurations, the I ++ and I data indeed lie 
together for the same magnitude of In , and the suppres- 
sion of switching threshold for this pair is higher than for 
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FIG. 3: Results related to the second moment of shot noise, 
(a) Effective temperature extracted from measured escape his- 
tograms using Eq. (2). At each noise current T* has been 
evaluated at various values of escape probability in the range 
0.1 - 0.9. The solid line is obtained using Eq. (JTJ with F = 1 
and Q — 2.5 determined from the bias current at which the 
underdamped phase diffusion sets in. This is shown by the 
arrow in (b), where width of the histogram has been plotted 
against In- The dashed line is from the effective temperature 
model, ignoring phase diffusion. 



P and I '". The latter two are again grouped together, 
as depicted in Fig. (b) . This shift is the manifestation 
of the skewness of the current distribution, where the 
sign of the third moment is positive. The values of Iq 
and C of the detector in Table I are obtained from fits 
to MQT escape [H; Fig. [2(c) shows the rates. The ob- 
tained value of Iq agrees with the Ambegaokar-Baratoff 
prediction [2l[ , and C with that based on the size of the 
junction. 

In Fig. [3] we have analyzed the escape histograms of 
Fig. Abased on Eq. © to extract effective temperature 
at various levels of In- We obtain in (a) a set of val- 
ues of T*, spread within about ±2% around the mean, 
for each In at different values of bias current I along a 
histogram. For clarity, we show data of one pulse con- 
figuration (7 H ) only, because the third order effects are 
weak: the variation in T* between different configura- 
tions is only about 1% at maximum. Figure [3] (b) shows 
the measured width, defined as AI between P — 0.9 and 
P = 0.1 of the same escape histograms. We see an in- 
crease up to about 7n = 2.6 fiA, in fair agreement with 
the thermal activation model (dashed line). Above this 
value of 7n the width starts to decrease. We attribute the 
maximum to the crossover from escape dynamics to UPD 



17| , now induced by the increase of effective temperature 
and subsequent decrease of the switching threshold cur- 
rent. From the value of I/Ic — 0.52 for the P = 0.5 
switching threshold at this noise current, we obtain the 
quality factor Q = ^f- ~ 2.5. (Similar procedure gave 
Q = 2.5 also for Sample SIS.) Based on this value of Q 
and the values of 7c and C from independent measure- 
ments, we are then able to draw the theoretical line based 
on Eq. Q in Fig. [3(a). 

Figure f?] shows the influence of the odd moments on 
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abatically fluctuating bias currents. A similar problem 



FIG. 4: Difference in switching currents of the threshold de- 
tector under reversal of the relative polarities of bias and noise 
currents. Data of Sample NIS are shown by filled dots, of 
Sample SIS by triangles and data of the reference sample, 
Sample REF, by open dots. The solid lines are the result 
of the theoretical model for Samples NIS and SIS. At cur- 
rents 7n < 0.2 jj,A the noise source of Sample SIS is in su- 
perconducting state, and the signal due to the third moment 
vanishes. 



has been analyzed theoretically in 11| using an effective 
Fokker-Planck method and very recently in [3] by gen- 
eralized stochastic path integral methods. The results of 
our model coincide with ll| up to numerical prefactor 
of the order of unity in the limit of large Q, and with 
[l3| excluding the circuit corrections that are weak in 
our case. Finally, although the data presented here can 
be interpreted quantitatively, it remains a challenge to 
engineer the radio-frequency environment on the chip to 
the degree that no uncertainty in the relevant bandwidth 
would remain. However, it is already possible to envision 
experiments where different mesoscopic noise sources use 
a common J J detector of third order fluctuations, and 
where a tunnel junction serves as a reference scatterer. 

We thank J. Ankerhold, H. Grabert, T. Ojanen, H. 
Pothier and E. Sukhorukov for useful discussions and 
Academy of Finland for financial support. A.V.T. thanks 
Finnish Academy of Science and Letters for a scholarship. 



escape threshold of all the three samples. The quantity 
1+ - I~ has been extracted as + I~ + ) - (I ++ + 

I )], all taken at P — 0.5. We see that the signal due 
to the third moment, plotted against Jn again, of Sample 
NIS (filled dots) has a maximum value of about 0.6% of 
Ic- This is to be contrasted to the essentially vanishing 
signal of the reference sample (open circles). The solid 
line following closely the data of Sample NIS is the re- 
sult of the theoretical model according to Eqs. ([3|) and 
(|4"|) . The non-trivial parameters needed to construct this 
curve, dl/dP and T* , have been extracted from experi- 
ment. Note that values of F and Q are not needed here 
since we measure T* directly. In addition, we have set 
F3 — 1, the value for Poissonian noise of a tunnel junc- 
tion, and Aw = w p , which is a very natural choice for 
our model including adiabatic frequencies only. We also 
show similar data of Sample SIS (triangles). Here, the 
bandwidth was taken to be 0.8u> p for the theory line. It 
may seem surprising at the first view that the results fol- 
low the theoretical prediction also in the regime where 
the detector is in the UPD regime (In > 2.6 /iA for Sam- 
ple NIS). Yet this is quite natural, because the frequency 
of the phase diffusion events is < 100 MHz, whereas the 
relevant bandwidth of noise affecting the phase particle 
extends up to ~ 50 GHz. Therefore it suffices to study 
escape from one well only. 

We have observed the third moment of shot noise up 
to frequencies of ~ 50 GHz, and interpreted this obser- 
vation by assuming that the third order effects manifest 
themselves in our JJ system as an average thermal es- 
cape rate in a potential which varies according to adi- 
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